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I. INTRODUCTION 

In the recent paper^, the scaling theory of electron lo- 
calization is discussed. It is argued that the standard 
interpretation of numerical data based on the finite size 
scaling analysis 2 4 is not correct. For the quasi-one di- 
mensional Anderson model, new formulation of the scal- 
ing, based on the analytical self-consistent theory, is pre- 
sented. The theory gives for the three dimensional (3D) 
Anderson model the critical exponent v = 1, in agree- 
ment with original self-consistent theory of Anderson 
localizatiorP. New scaling relations have been proposed 
for higher dimension d > 4. 

In this comment we show that the theory^ is not con- 
sistent with present numerical data for the 3D and 5D 
Anderson model. 

We consider Anderson modeP with diagonal disorder 
W defined on the quasi-one dimensional system of the 
size 



L 



d-l 



x L, 



(1) 



{d is the dimension of the model) and calculate the small- 
est Lyapunov exponent z\{W,L). The later is related to 
the localization length 



2L 



zi = 



;id 



(2) 



and determines the exponential decrease of the wave 
function, |*| 2 ~ exp[-z 1 L z / L}. 4 For the 3D model, 
L z = 2L/e 2 is sufficient to achieve the relative numer- 
ical accuracy The size L varies from L = 8 to L = 34 
for d = 3 and is L < 8 for d = 5. 



II. 3D SYSTEM 

Suslov's theory predicts^ that in the vicinity of the 
critical point (r = W — W c Cl) the localization length 
follows the scaling behavior 
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= y* + At(L + L ) 



(3) 



with a new additional length scale Lq not considered in 
the standard scaling analysis, (y* is the size-independent 



critical value). This prediction is in variance with the 
standard scaling formula 



~i 



2L 



Zlc + 



(4) 



used in the finite size scaling analysis of numerical 
data.™ 

To support the result ([3]), Suslov used numerical data 
for parameter Z\ published in RefP and found that L rj 
5 (left Fig. 6 irP). We show in Fig. [ljthe same Figure 
with additional data for 24 < L < 34. Power fit Z\(V) = 
a + bL a calculated for W — 16 and W = 17 supports the 
validity of the the relation Q . 

Before testing the validity of Eq. Q we have to notice 
the relation ^ between the localization length expressed 
in Eq. ([3| and the parameter z\ shown in Fig. [I] We fit 
our data for z\ to the function 



c = 



1 



a + a\L 



(5) 



shown by dotted lines in Fig. [T] Comparing with Eq 



3} and using y* = z^ = 3.48" 1 (Fig. V2j we obtain 
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FIG. 1. The 3D Anderson model: The parameter z\{Li) for 
various disorder Solid lines are power fits for W — 16 and 
W = 17. Contrary tcP, fits are not linear in L. Note that Z\ 
decreases for W = 16.5 and increases for W = 16.6. There- 
fore, we expect that 16.5 < W c < 16.6. Scaling analysis 
gives W c ~ 16.55. Dotted lines are fits (5} with a = 0.302 
and ai = 0.0017 (W = 16) and a = 0.267, ai = -0.00108 
(W = 17). 
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L = 20 3.48 + 0.97 L + 0.07 L 2 

L = 26 3.48 + 1.17 L + 0.19 L 2 

L = 30 3.49 + 1.29 L + 0.17 L 2 
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L = 34 3.48 + 1.36 L + 0.20 L 2 
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FIG. 2. Quadratic fit |6| of z\(yV — W c ) for four values of the 
size L. The cross section determines z\ c ~ 3.48. 
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FIG. 3. The L-dependence of the slope s(L) ~ L 1/u . The 
critical exponent v = 1.566. Dashed line shows the linear 
L-dependence, predicted by Eq. pj. 



Lq ps 8.6 from W — 16 data, but significantly different 
value L « 17 for W = 17. 

Although the power fit Q is clearly better than the 
fit ([5]) , Fig. [I] shows that the estimation of true scaling 
behavior might be difficult since various analytical func- 
tions seem to fit numerical data with sufficient accuracy. 
In the present case, the problem lies in the non-zero crit- 
ical value Z\ c . To avoid the ambiguity in the choice of 
the fitting function, we have to extract the critical value 
from numerical data^. When data for Z\ are plotted as 
a function of the disorder (Fig. [2]), we can fit them by 
quadratic polynomial 

Zl (W,L) = z lc + Ts(L)+T 2 t(L) (6) 

and calculate the L-dependence of the slope s(L). From 
Eq. pj we see that s(L) should be a linear function of L, 
while Eq. Q predicts power-law behavior s(L) ~ L x l v . 
Figure [3] shows s(L) as a function of L. The fit confirms 
the power-law dependence s(L) ~ L 1 /" with critical ex- 
ponent v m 1.56, as obtained by other method^. 

III. 5D MODEL 

For higher dimension, the following size dependence of 
the localization length at the critical point (r = 0) was 
derived 

t / t \ (d-4)/3 

¥ - (;) ■ < 7 » 

In particular, for d = 5 Eq. Q gives 

z 1 (r = 0)^L- 1 /3 ( 8 ) 

which means that the critical value of Z\ is not size inde- 
pendent but decreases to zero when L — > 00. Since the 
localization length is finite in for t > 0, the r-dependence 
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FIG. 4. The parameter zi(W) for various system size. For 
W < 57.5, z\ decreases when L increases. There is no indica- 
tion for the critical behavior described by Eq. Q 

of Z\(L, r) for fixed L must exhibit an infinite disconti- 
nuity at r = 0: 

Zl(r) ~\L r>0 (9) 

We test numerically the size and disorder dependence of 
z\ . We show in Fig. [4] and [5] the disorder dependence 
of z\ for fixed L. Our data in Fig. [4] do not indicate 
any discontinuity in the L dependence. Contrary, z\ is 
smooth analytical function of both parameters, W and 
L. 

For smaller disorder, z\ is always decreasing function 
of L. This is typical for the metallic regime. However, 
Z\ does not depend on the size L when W = 57.5. This 
is consistent with the scaling equation Q. Insulating 
regime, where Z\ increases with the size L is observed 
only when W > 57.5 (Fig. [5). 

Note that z\ ~ 7 for disorder W ~ 57.5. Therefore the 
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FIG. 5. The 5D Anderson model: the parameter Z\ as a 
function of disorder W for L = 4, 5, 6, 7 and L = 8. Data 
indicate that z\ does not depend on the size L when W ~ 
57.5. This value is considered as a critical disorder W c in 
the "standard" finite size scaling theory. Solid lines are fits 
Zi(L) = Zi + s(L)(W — Wc). Inset shows the //-dependence 
of the slope s(L) ~ L 10413 . The original figure was published 
irP but new data for L = 8 were added. 



localization length, 

6d = —L (10) 

is much smaller than the size of the system and we do 
not expect that finite size effects play significant role al- 
though the size L is much smaller than in 3D system. 

Scaling analysis, similar to that for the 3D model en- 
ables us to find the critical exponent, rj 0.96. 

IV. CONCLUSION 

We showed that numerical data for the parameter z\ 
do not agree with the predictions of the theory^. Both 
Z\ and the localization length are analytical continuous 
functions of the disorder W and the size of the system L. 

For the 3D system, we presented additional numeri- 
cal data for larger system size L < 24 up to L = 34. 
These new data confirm previous estimation of the criti- 
cal exponent v — 1.56!^ It is worth to mention that the 
same value of the critical exponent was obtained already 
20 years ago with the use of numerical data for L < 12 
only^l. w e also note that the same value of the critical 
exponent was obtained from numerical analysis of other 
physical quantities: mean conductance, conductance dis- 
tribution, inverse participatio n rati c^ and also for critical 
points outside the band center^^l. This value of critical 
exponent was recently verified experimentally^ and cal- 
culated analytically-^. 
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